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Abstract 

The scalar dark matter candidate in a prototypical theory space little Higgs model is 
investigated. We review all details of the model pertinent to a relic density calculation. We 
perform a thermal relic density calculation including couplings to the gauge and Higgs sectors 
of the model. We find two regions of parameter space that give acceptable dark matter 
abundances. The first region has a dark matter candidate with a mass 0(100 GeV), the second 
region has a candidate with a mass greater than O(500 GeV). The dark matter candidate in 
either region is an admixture of an SU(2) triplet and an SU (2) singlet, thereby constituting a 
possible WIMP (weakly interacting massive particle). 



1 Introduction 



The "little Higgs" models provide a new mechanism to stabilize the weak scale wherein the Stan- 
dard Model Higgs is a pseudo Goldstone boson |TJ E3 H|- The mass of the Higgs is protected by 
approximate non-linear global symmetries. Subsets of these global symmetries are broken by cou- 
plings in the theory; however, this breaking occurs in such a way that any single coupling preserves 
sufficient symmetry to keep the Higgs massless. This ensures that there are no one loop quadratic 
divergences to the Higgs mass and allows a two loop separation between the weak scale and new 
physics. The properties of these models allow for some distinct phenomenological signatures El - 

The dark matter problem is one that has long been optimistically tied to the theory of weak 
interactions. At present time, it is not certain exactly what constitutes the dominant portion of the 
matter density in the universe, but its non-baryonic nature is generally agreed upon. It happens 
that a stable neutral particle with Standard Model weak couplings and a weak scale mass gives 
roughly the correct amount of thermal dark matter. Such an intriguing coincidence should not be 
overlooked in any theory attempting to explain the breaking of electroweak symmetry. Since little 
Higgs (LH) models offer a mechanism to stabilize the weak scale, it is natural to ask if LH models 
have anything to say about the dark matter problem. 

In this paper we investigate the possibility that LH models may explain the observed abundance 
of dark matter. We first review some general features of theory space LH models that are related 
to a dark matter analysis. Then we specialize to the case of the 50(5) theory space LH model 
contained in Next we perform a relic abundance calculation including couplings of the dark 
matter candidate to the gauge and Higgs sectors of the model. We find that a specific LH model 
allows proper dark matter abundances in two distinct interesting regions. First, the correct relic 
abundance is achieved when the dark matter candidate, Ni, has mass (9(100 GeV). A mass of 
(9(100 GeV) is allowed for either a weakly coupled or a 'super- weakly' coupled iVi. By super-weak 
coupling we mean that it does not couple significantly to Standard Model (SM) weak vector bosons, 
it only couples to the TeV scale vector bosons. A second preferred regime exists where the dark 
matter candidate is heavy, with a mass greater than (9(500 GeV). The preferred couplings of Ni 
in this regime vary from weak (at of (9(1 TeV)) to 'super-weak' (at rriN 1 of (9(500 GeV)). We 
discuss possible additional effects on the relic density calculation (beyond couplings to the gauge 
and Higgs sectors) and find that they will usually only strengthen our conclusions. 

2 Dark Matter Candidates in Little Higgs Models 

In terms of dark matter, LH models come in several varieties. Some LH models, referred to as 
'theory space' LH models [D 121 IS], contain new possibly stable neutral particles. Many other LH 
models, especially non-theory space models jlj , have no potential dark matter candidates associated 
with the little Higgs mechanism. For those theories that do contain the possibility of new stable 
neutral particles, one can place constraints on the LH model by requiring that this particle supply 
all of the necessary dark matter. In the very least, one must require that this particle not provide 
more than the required dark matter. If it provides less, then to explain the remaining dark matter 
there must be other stable particles from the ultraviolet completion to the LH model, potentially 
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similar to particles in 

At this stage it is useful to state what specific characteristics of LH models allow for a possible 
dark matter candidate and to compare these with models of broken super symmetry. In models of 
broken supersymmetry, it is usually R-parity that keeps the lightest supersymmetric particle (LSP) 
stable. One way to define the charge of a given particle under R-parity is: R = (— i) 2 i+ 3B + L . 
Here j is the spin of the particle, B is the baryon number of the particle and L is the lepton 
number of the particle. This results in every ordinary particle being even under R-parity and every 
supersymmetric partner being odd under R-parity. If R-parity is a conserved quantity, the lightest 
R-odd particle must be stable. This stability is one of the crucial ingredients that allow the LSP in 
models of broken supersymmetry to be a possible dark matter candidate. Since the present number 
density of a thermal dark matter candidate is determined by the evolution of its abundance on a 
cosmological timescale, any sizable violation of R-parity will result in no observable present-day 
relic density. The other necessary ingredient for dark matter is the correct combination of mass 
and annihilation cross section. It is an intriguing, and much-studied, coincidence that a stable 
particle with a mass of order the weak scale and annhilation rate dominated by weak processes 
gives a present thermal relic abundance in the right range to explain the observed dark matter. In 
models of broken supersymmetry, the lightest neutralino and sneutrino have just these additional 
characteristics, on top of one of them frequently being the LSP. 

What characteristics of LH models might allow for the right dark matter abundance? Are the 
requirements of stability, weak-scale masses, and weak coupling accurate and sufficient to provide 
the correct relic density? We will find the necessity of weak-scale masses and weak coupling to 
be model-dependent. However, let us repeat that the issue of stability is generally not a model 
dependent question. Since cosmological timescales govern the annihilation (and possible decay) of 
a dark matter candidate, any such candidate with a lifetime shorter than the age of the universe will 
have completely decayed by today, resulting in no relic abundance. Fortunately some LH models 
do contain symmetries that can protect the stability properties of some constituent fields. 

Theory space LH models allow for a geometrical description of the group and field content. Such 
a geometrical description can also be used to illustrate some of the symmetries of the models. For 
instance, many theory space LH models contain a Z 4 symmetry that can interchange the non-linear 
sigma model fields among themselves. This symmetry is visualized in [2 j as rotations of 90° on a 
square torus, though the symmetry exists in many theory space LH models without such a simple 
visualization. For instance, this symmetry is contained in the model 7\ we study in greater detail 
later on in this paper. Because the lightest Higgs boson is contained in these non-linear sigma 
model fields, the Z 4 symmetry is broken once electroweak symmetry is broken by the development 
of a Higgs vacuum expectation value (vev). In the model we study in detail, the existence of the 
Higgs vev breaks this symmetry down to Z2. It is the existence of this symmetry that allows theory 
space LH models to contain interesting dark matter candidates. 

Given the freedom contained in many LH models, it is possible that the Z 4 and consequent Z 2 
symmetries be approximate instead of exact. This is disastrous for the dark matter candidates we 
investigate, as has been explained above. Therefore, our results apply only to models in which these 
symmetries can be made exact. From here on we make the explicit assumption that at least the Z2 
symmetry is exact. As a result of our assumption, the lightest particle charged under this Z2 will be 
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stable. In theory space LH models, the Standard Model fields, including the Higgs, are even under 
this Z 2 . In addition, theory space LH models always contain scalars and rj, which are both odd 
under the Z 2 . Thus, the lightest of these two will be stable and will constitute the model's only hope 
for a thermal dark matter candidate. Interestingly enough, and r\ are uncolored and also transform 
under the electroweak SU(2)w x U(l)y as 3o and lo, respectively. So only the can be considered 
a WIMP. Regardless, there are important similarities between particle characteristics that exist in 
theory space LH models and characteristics that allow the LSP in broken supersymmetric theories 
to provide roughly the correct amount of dark matter. 

After electroweak symmetry breaking (EWSB), and T] decompose into 

0^(0°,0±) T] ^ T] o (1) 

Typically the 0° and f]° mix after EWSB due to interactions of the form rjh'(j)h. We label the 
mass eigenstates Aq and N[. Electrically charged particles cannot serve as dark matter so one must 
require m^± > m^ 1 . The and rj can either have the masses generated completely radiatively or 
have a tree-level mass from an U Q plaquette" Schematically, the mass matrix is of the form: 

^Mass = -krnlM ) 2 + ml Srff-lm* 0V)-mL|0 ± | 2 (2) 



with 



m H = m G + m l + m E + m l 
2 2,2,2 

m m = m s + m E + m n 

m ln = m l (3) 

where tuq, m|, m% and, ttlq are the masses from gauge interactions, scalar interactions, electroweak 
symmetry breaking, and possibly a tree-level mass from an Q plaquette. Because rj is uncharged 
under the Standard Model gauge interactions, it does not receive a mass from this sector of the 
theory. The gauge contribution is typically the largest radiative contribution to the masses because 



2 3g 2 / A 2 \ 2 A 2] AM 2 1 x 2 

lG * — 2 m wl \o g ^-2- j m s * — rVog ( ) m B * -Xv , 



with mw' ^ 1-8 TeV and ~ 1.5 TeV 1 . Here 0' are the heavy scalars needed to cancel the 
quadratic divergences of the Higgs quartic coupling and W are the heavy gauge bosons needed to 
cancel the divergences coming from the Standard Model gauge bosons. The parameter A determines 
the weight of the Higgs potential to be defined in Eqn. [T7| As usual, v is the magnitude of the two 
Higgs vevs, v 2 = v \ + v\ for the two Higgs doublet model we study in detail in the next section. 
The cut-off of the theory is A ~ 4irf ~ 10 TeV and g is the gauge coupling related to the TeV scale 
vector bosons. The mass matrix for 0° and rf can be diagonalized by an orthogonal transformation 



1 The constraints on the masses of these particles come from precision electroweak constraints and limits on the 
breaking scale / 0U2EI]. 
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with a mixing angle, i?^. If contributions from f2 plaquettes are ignored, the lightest eigenvalue of 
the mass matrix is of order 100 GeV. Including Q plaquettes can lift the lightest eigenvalue up into 
the TeV range. 

Depending on the full model, there can be a symmetry that interchanges <f)° <-> rf in the 
absence of gauge interactions. For the SO (5) model we study in the next section, this symmetry 
just exchanges the SU{2)i and SU{2) r subgroups of the 5*0(5). This symmetry guarantees that the 
scalar (nig) and EWSB (m 2 E ) contributions between <p° and 77° exactly cancel in the calculation of 
this mixing angle, now given by 



2m 2 

tan 20^ = —^. (5) 



This angle becomes small when w? G grows large. In this limit the dark matter candidate becomes 
predominantly ?7°-like and does not have gauge interactions with the electroweak vector bosons. 
However, rf does have gauge interactions with the additional vector bosons of TeV-scale mass. It is 
for this reason we call this limit the "super- weakly" interacting limit. It is not uncommon to have 
tuq large enough to yield cos 2 ^^ ~ 0.95. 



2.1 The Minimal Moose 

The smallest theory space LH model is called the "Minimal Moose" |J]. The total global symmetry 
structure of the theory is (Gl x GrY- This global symmetry structure is the essential feature that 
protects the Higgs from one loop divergences. The theory contains four non-linear sigma model 
(nlerm) fields Xi, i — 1, ...,4. The presence of each nlamfield breaks a global Gl x Gr symmetry 
down to the diagonal group Go- Under global symmetry transformations and gn, the nkrm fields 
transform as 

X t -> g Li Xig Ri . (6) 

In some models there is also an exact Z4 discrete symmetry that interchanges the nkrm fields. It is 
the Z4 symmetry that will result in a stable dark matter candidate. 

Only a subgroup of the entire global symmetry group is gauged. Inside G R , 577(2) x U(l) is 
gauged while inside G L , a group G is gauged which also contains an 577(2) x U(l) subgroup. The 
four nkrmfields break this extended gauge symmetry down to the diagonal SU(2) W x U(l)y. The 
nlam fields transform under gauge transformations as bi-fundamentals: 

Xi -> gGXigl u{2)xU{iy (7) 

In terms of the vocabulary of 'moose' diagrams, this setup is a 'two site - four link' model with each 
'site' being a gauge group (G and 577(2) x £7(1), respectively ), and each 'link' being a nlcrm field. 
Since each link is charged under the gauge groups living at both sites, a geometrical picture can be 
drawn of the model as two dots connected by four lines. This model is closely related to the 2x2 
toroidal moose of [T], which is a four site model and is also easy to visualize. One views the 2x2 
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toroidal moose as four dots, the corners of a square, connected by link fields to make the structure 
of a toroid. The minimal moose is related to the 2x2 toroidal moose through orbifolding by a 
translation along each of the diagonals of the 2x2 toroidal moose. 

The gauge symmetries explicitly break the global symmetries. The breaking of the G\ global 
symmetry is accomplished by only the G gauge transformations. Similarly, the breaking of the 
G A R global symmetry comes only from the 577(2) x £7(1) gauge transformations. Since the gauge 
structure of the model is known, one can write down the kinetic terms in the Lagrangian: 



where / is the 'pion decay constant' of the nkrm. The are the covariant derivatives, to be 
defined more explicitly below. 

Up until this point we have kept the gauge group G unspecified. We have done this because the 
minimal moose is, in some sense, modular. Many general properties are determined irrespective of 
which group is chosen, as long as G contains 577(2) x £7(1). However, each group offers different 
predictions for TeV-scale physics. As a result, some choices lead to tension with measurements of 
precision electroweak observables. For instance, in the original minimal moose G was chosen 
to be SU(3). The theory resulting from this choice can be significantly constrained by precision 
electroweak physics [TJ]. In [JJ it was shown that SO (5) (or equivalently 5p(4)) has a custodial 
SU{2) symmetry that allows a simple limit where precision electroweak constraints are rather easily 
satisfied. There are several other groups that would work equally well such as S77(4) or SO (7), 
but we choose here to only pursue the SO (5) model described in |7j. This group has the added 
bonus that there is a symmetry in the gaugeless limit that interchanges <f)° and rj so that Eq. |3] 
is valid. Another important motivation for SO (5) is that the approximate Z 4 symmetry that is 
phenomenologically necessary to prevent unacceptably large triplet vevs can, and will, be lifted to 
an exact symmetry and will ensure the stability of the dark matter candidate. 

Even though the SO (5) model is relatively unconstrained by precision electroweak physics, 
some constraints from such considerations do exist. The most important comes from the breaking 
of custodial 577(2) by electroweak symmetry breaking vevs. This breaking of custodial 577(2) places 
a limit of 700 GeV on the nlam breaking scale, /. The custodial 577(2) violation is proportional 
to sin 2/3, where tan/5 = v%/vi is the ratio of the vevs of the two Higgs present in the model (see 
below). To suppress this effect, it is preferred for (3 to be relatively small. In the Z4 symmetric 
limit that we are pursuing to find dark matter, there are no issues of triplet vevs because there are 
no trilinear couplings between the Higgs doublets and the SU(2)w triplets. The effects upon the 
oblique S and T parameters from the new Z4 symmetric top sector would need to be calculated in 
this model. However, it should be possible to make the effects small without significantly altering 
the limits. In addition to the mildness of precision electroweak constraints, another benefit of this 
model is that it has a relatively minimal set of particles. 

Let us explicitly summarize some details about the symmetries and symmetry breaking in 
the model we are investigating. The SO (5) minimal moose defined in [JJ has the expected site 
group structure G x G' = 50(5) x [577(2) x £7(1)]. The [50(5)] 8 global symmetry is broken by the 
50(5) x [SU (2) x £7(1)] gauge interactions. The aforementioned custodial SU (2) is approximate and 




(8) 
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comes about from the gauging of 5£7(2) r in addition to the secondary SU (2)i, both from the 50(5). 
Thus, there exists a complete set of 5*0(5) vector bosons with SU(2)i x 5£7(2) r representations 
W l ~ (3,, l r ), W r ~ (l b 3 r ), V ~ (2,, 2 r ). Due to the embedding of U(1) Y inside the full 50(5), 
the W r bosons split into W r,± and W r ' 3 . W r ' 3 is the field responsible for cancelling the quadratic 
divergence of the £7(l)y gauge boson. If it is assumed that the 50(5) gauge coupling is large, the 
Standard Model W and B are mainly composed of the [577(2) x £7(1)] gauge bosons. The heavier 
orthogonal combinations, W and B' are thus mainly composed of the 50(5) gauge bosons. The 
nlam fields X { break 50(5) x [SU(2) x £7(1)] down to the Standard Model SU(2) W x £7(l) y , thereby 
giving W and B' masses in the TeV range. At scale of the W mass, the entire Z4 symmetry is still 
intact. As the energy scale is lowered, the standard Higgs mechanism breaks electroweak symmetry 
and the Higgs vevs break Z 4 — > Z 2 . However, as previously mentioned, this Z 2 is sufficient to 
preserve the stability of the lightest eigenstate of the (77, </>) system. 

In the case of 50(5) the covariant derivative acting on the nkrm fields is 

= d^Xi - ig 5 X t T^w!™ n] + 1 (g 2 T l ^ a + 9l T r > 3 W r /) X,. (9) 

In this definition, are the gauge bosons of the group 50(5), W l ^ a are the 5£7(2) gauge bosons 

and W^' 3 is the £7(1) gauge boson. Expanding the nkrm fields around Aj = 1, one finds that the 
kinetic term contains mass terms for the TeV-scale vector bosons with 

a % 2 / 2 2 8g'*p 

mw ' = sim^' mB> = s^2e>- (10) 

The mixing angles 9 and 9' are defined as 

tan# = ^, taji9' = —. (11) 

95 95 

At low energy the minimal moose has the same physics as any toroidal theory space jT| 
At tree level, two orthogonal combinations of pseudo-Goldstone boson multiplets are massless. In 
terms of the linearized modes A, = exp(ixj//), we take the massless combinations to be x\ — £3 
and X2 — X4. Continuing to follow the conventions of [3], these massless fields can be parametrized 

as 



X l = Xl = e 2i{x+y)/f (12) 

X 2 = X\ = e 2i{x - y)/f . (13) 

So in terms of the massless modes, there are two nlcrm fields, A = exp(ix/ f) and Y = exp(iy / /) , 

called the little Higgses, that transform as adjoints under the diagonal global 50(5) symmetry. The 
low-energy Lagrangian in terms of these fields is: 

f 2 f 2 

Cle = yTrl^Al 2 + ^-Tr\D^Y\ 2 + A/ 4 TrAFA t F t . (14) 
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The covariant derivative, D, acting on X and Y, only contains the Standard Model SU(2)w x U(l)y 
gauge fields. X and Y also have interactions with the TeV-scale W and B' vector bosons, but these 
have not been included in Eq. El For dark matter calculations, the most important interactions 
are given by Eq. El and 

Cw> int = ig cot 26 (TtW'^D^X + TrW'^D^ + h.c. + • • • (15) 
£b> int = ig' cot 26' (TtB'^D^X + TrB'^D^ + h.c. + • • • , (16) 

where interactions involving more than one heavy vector boson have been suppressed. The gauge 
mixing angles 6 and 6' have been given in Eq. ^2 Bounds from non-oblique precision electroweak 
corrections require both 6 and 6' to be small FEU [7J EJ E] • For simplicity we take the couplings 
of X and Y to W and B' to be identical to the couplings to the Standard Model W and B. 

The potential in Eq. El when expanded in terms of the linearized modes to leading order, is: 

V(x,y) = \Tr[x,y] 2 . (17) 

The linearized fields, x and y, are adjoints under the unbroken global Gd- Under the electroweak 
SU{2) W x U(1)y C Go-, the linearized fields decompose into several representations. These rep- 
resentations include and 77, which are in the 3o and lo representations of SU(2)w x U(l)y, 
respectively. The decomposition also contains h, which is in the 2i representation and is identi- 
fied as the Standard Model Higgs. There can be additional matter fields as well, but those are 
unimportant for our present discussion. 

The discrete Z 4 symmetry of the low energy theory in Eq. ITU takes 

X -> Y, Y -»■ Xl (18) 

In j2] this symmetry is understood as 90° rotations on a square torus. It is possible to lift this 
symmetry to the entire high energy theory and also to make the symmetry exact, as is done in 
In many models there will be a transformation that takes X — > QXQ and leaves the 77 and fields 
invariant but rotates the Higgs fields by e l7T . For instance, such a transformation can take the 
form of Q = diag(— 1, — 1, 1) for G = SU(3) [2] or Q = diag(— 1, — 1, — 1, — 1, 1) for our choice of 
G = SO (5) [7j. For the model of this is just the Z 2 symmetry from the SU(2)i subgroup of 
SO(5), different from the Z4 symmetry of Eq. El After electroweak symmetry breaking the discrete 
Z 4 x Z 2 structure breaks leaving only an exact Z 2 symmetry: 

X -> QX^Q Y -> QY^Q. (19) 

Under this symmetry, the Standard Model fermions and Higgses are even, while <fi — ► —<p and 
r) — > —r]. This symmetry will leave the lightest triplet or singlet scalar stable and can result in dark 
matter. It is important to note that to have a viable dark matter candidate it is necessary both 
that the Z 2 be exact and that the dark matter candidate be the lightest field transforming under 
this Z 2 . Generally the stable scalars in X do not mix significantly with the stable scalars in Y. This 
is because mixing is only generated at loop level through higher order interactions with the Higgs 
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and the top quark. Each little Higgs, X and Y, contains a neutral r\ and a neutral 0, so there are 
four electrically neutral mass eigenstates that come from r) Xjy and 4> XtV . The interactions of the X 
and Y fields are the same at tree level, so the minimal moose has two almost-identical copies of the 
neutral r\ and fields. Thus we typically have two nearly degenerate scalars N\ and N2 and the two 
heavier admixtures N[ and N' 2 . Since N\ and N2 are nearly degenerate, the evolution of their relic 
abundances are intertwined. Even though the heavier of the two (assumed to be N 2 for simplicity) 
will decay into the lighter one quickly (compared to the time scale involved in relic abundance 
evolution), both play important roles in the correct determination of the final dark matter density. 
When we perform the relic density analysis we take this complication into account. 

It is also possible to deform the minimal moose by adding Q plaquettes. At low energies these 
deformations result in terms being added to Eq. El of the form 

£n Piaq = nfTiVLXVLX + KfTrOYttY + h.c. . (20) 

These operators give tree-level masses for the scalars rj and while leaving the Higgs massless. 
They can be made to be symmetric under the full Z4 x Z2 high-energy discrete global symmetry. 
Furthermore, these deformations allow us to lift the masses for the 0s and rjs from the weak scale 
to the TeV scale without affecting naturalness. 

It is important to remember that the mass matrix of the (77, 0) system is dominated by contri- 
butions from gauge interactions, and this frequently results in a lightest eigenstate that is mainly 
//-like. Due to this strong dominance in determining the mass eigenstates, the gauge interactions 
also have a dominant influence on the relic abundance, as will be clearly shown in the next section. 
For concreteness, we explicitly give the part of the covariant derivative containing 0° : 

U,<p x>y D ^_ ^ %gW2A + J (21) 

The charged and neutral fields also have interactions with the heavy vector bosons through 
interactions derived from Eq. 

Cw mt = ig cot 26 (TrWfolDffa + TrW^j, + h.c. + • • • (22) 

Cb> mt = ig' cot 26' (TtB'^Id^ + TrP^D"^ + h.c. + • • • . (23) 

We have left the in its SU{2) triplet form in Eqs. 1221 an d 1221 but decomposition into 1 * 1 
and 0° follows using the standard operations. One can easily find the interactions involving 0° by 
realizing that the neutral field comes only from the a = 3 part of = a T a where T a are the SU{2) 
generators defined in the Appendix of Ref. [7j. 

The rf has no direct gauge interactions except for terms like g' 2 r] r] B' B' 1 * . These terms do not 
contribute to an extent that suppresses the relic density into the cosmologically preferred range. 
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After the gauge interactions, the most important interactions come from the quartic potential. 
To determine the quartic potential, one must define the vevs for the linearized modes x and y which 
contain the Higgses. The directions of these vevs are: 

M - VCOS l 3cOB 2 T vO V COS ^ Sil1 2 T v3 ( X _ VSill f 3cOS 2 T vO , V _ Sil1 g Sin 2 jn, 3 ^ 

where we are again using the generators in the Appendix of Ref. |7j. The quartic potential cannot 
stabilize EWSB when £ = 0, tt because this is precisely where [x,y] =0. In the angle £ was 
taken to be | to simplify the analysis, but it is possible that a symmetry can enforce this. In this 
paper we use this limit as a simplifying assumption. 

Below we give the interaction potential for the (rj, <p) system that we use in part to compute 
the annihilation rate in the next section. This interaction potential can be derived directly from 
Eq. El The h is the usual lightest Higgs, H° is the heavy CP-even Higgs and A is the CP-odd 
Higgs. The h° — H° mixing angle, a, is typically close to (3 with this Higgs potential although our 
results are generally insensitive to this. Choosing a = (3 the potential is: 

V = y />°(cos 2 m - tf + sm 2 Ml ~ ill?) 
+ ^sin2/3i/ o ((0°-r / O) 2 + (0;-^) 2 ) 

+ < 25 ) 

As can be seen in the above equation, it is only interactions with the CP-odd Higgs that allow the N\ 
and N2 eigenstates to coannihilate 2 . Other than these interactions, the cosmological development 
of their relic densities is independent. 

On the surface, the physical Higgs sector of this model resembles that in the MSSM. However, 
one notable exception is the coupling of b quarks to the lightest Higgs. In the MSSM, the coupling 
of the lightest Higgs to b has a fixed tan /3 dependence. In the minimal moose, there is a freedom 
as to whether the lightest Higgs couples to b quarks proportional to tan j3 or proportional to cot f3. 
This can have a non-negligible effect on the relic density. In the analysis below, we choose the 
coupling to be proportional to tan /3, but we indicate how the allowed LH parameter space would 
change had we chosen the coupling to be proportional to cot f3. 



3 Calculation of Relic Abundances 

Now we will use the general properties of the possible LH dark matter candidate to calculate its 
relic abundance. We will include all of the interactions listed in the previous section. These would 
be expected to be the dominant contributions. Subdominant contributions might be expected from 

2 Here we use the term coannihilation to refer to interactions between two different particle species in the thermal 
bath, such as N1N2 — > X . When we speak of standard annihilation, we mean interactions between two particles of 
the same species, such as N\N\ — > X. 
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interactions with other heavy states in the theory. We will find that a proper relic density generally 
prefers Ni to have a larger mass than estimated in the literature on LH models, though the notable 
exception exists of a super-weak eigenstate. 

We start here by assuming that have some masses mjvi = tun 2 an d can annihilate, as 
mentioned before, through a mixture of weak and super-weak gauge bosons and also the Higgs 
bosons 3 . We have given these interactions explicitly in Eqs. El IZH 1221 l221 an d|2ID From this we 
can calculate a thermally averaged cross section (erf). The thermally averaged cross section is used 
in the evolution equation for the number density of Nx, the Boltzmann equation: 



^ = -3Hn - (av) [n 2 - (n eq ) 2 ] (26) 

The Boltzmann equation shows that universal expansion always dilutes the number density of 
a particle species, but also that annihilations only become important once the particle species has 
diverted from equilibrium with the rest of the universe. One would naturally expect this departure 
from thermal equilibrium to happen once the temperature of the universe has dropped below the 
mass of the particle, T < m^. Once the particle has decoupled from the surrounding thermal bath, 
the number density is still being diluted by Hubble expansion. This expansion eventually shuts off 
the annihilation once the expansion rate, set by H, dominates over the annihilation rate, set by 
T = n(av). After this event occurs, called 'freeze-out,' the number density of the particle species 
is essentially frozen except for further dilution from Hubble expansion. Thus, one of the most 
important calculations to be done to determine a relic density is the calculation of the temperature 
at which freeze-out occurs, Tp. This temperature is determined iteratively through its dimensionless 
inverse xp = m^/Tp: 

x F = \n( *M*9rn Pl m Nl (<rv) \ (2?) 



y/9* %F 

Here g is the number of degrees of freedom for the dark matter candidate, rripi is the Planck mass, 
and g* is the number of effective relativistic degrees of freedom at the time of freeze-out of the dark 
matter particle. For weakly interacting cold dark matter candidates, one normally finds xf ~ 20. 
Having determined (av) and x F , the most important formula we need is 

Qdmh 2 = — y/2 T, — 7" ' ( 28 ) 

gj m-piJ {xp) 

where J(xf) = f^(av) x~ 2 dx. 

In many supersymmetric theories, one finds rridm — 100 GeV. Recall that in supersymmetric 
theories, many of the annihilation channels of neutralinos are suppressed by a velocity factor of 
(1 — Am 2 dm /s) because the initial state particles are fermions. Here s is the squared center-of-mass 
energy. This suppression is even more severe in neutralinos because they are Majorana fermions, 

3 We remind the reader that we also include coannihilation effects since N2 is assumed to be degenerate in mass 
with N\. These effects are explained in more detail in Section IXT1 
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Figure 1: Relic Density Contours for LH Models. The preferred relic density regions are 
displayed for tan/3 = 0.4. The green region illustrates a preferred value of 0.094 < VL dm h 2 < 0.129 
and the red regions are regions where dark matter is over-produced, Vt dm h 2 > 0.129. 



but the effect exists in some annihilation channels of Dirac fermions as well. Because the dark 
matter candidates are often scalars in theory space LH models, this suppression factor is absent. 
Therefore, the annihilation is more efficient and higher masses are required to achieve the correct 
relic density than in models involving neutralino dark matter. This absence of velocity suppression 
is also important in other models with scalar dark matter candidates [TEJ H3 UH UH 1201 EH 1221 123] 
such as inelastic sneutrino dark matter [21]. The heavier required masses can clearly be seen in 
Figure H] where contours of constant relic density are plotted in the (cos 2 1?^, m Nl ) plane. Here 
cos 2 "d^ is the mixing parameter defined in the discussion leading up to Eq. It is the mixing 
parameter that determines the admixture of weak and super-weak strength with which N\ can 
annihilate. In this and following plots, the green (light shaded) regions are cosmologically preferred 
with 0.094 < Qdmh 2 < 0.129 [37]. The red (dark shaded) regions are cosmologically excluded, 
having a relic density Qdmh 2 > 0.129. In Figs. Hand [2] we have chosen tan/3 = 0.4, a moderate 
value that is still allowed by precision electroweak constraints [7J. We have also fixed the lightest 
Higgs mass to be m h o = 136 GeV, the heavy Higgs masses to be m H o = m H ± = rriA = 500 GeV, 
and the heavy super- weak gauge boson masses to be mw> = ttib 1 = 1800 GeV. Furthermore, here 
we take the Higgs boson to couple to the b quark proportional to tan (3. As mentioned previously, 
this is an arbitrary choice. We explain below what happens if this choice is instead changed to 
cot (3. 

In Figure Q we can see that requiring the LH model to supply all required dark matter allows 
mjVi ~ 100 GeV for any value of cos 2 fi^, thus signalling a robust weak-scale dark matter candidate. 
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However, it is important to note that cos 2 ^^ > 0.8 is required to have m Nl not be constrained to 
lie on m h o/2 ± 5 or mw, where 5 is roughly 3 GeV. This behavior will be explained in the next 
paragraph. Between mw and m h o, a preferred region exists with the aforementioned cos 2 > 0.8. 
Above rriho, annihilation to an hh final state again reduces the relic density beyond relevancy. Once 
above m h o, one must go beyond mjVi ~ 500 GeV in order to find cosmo logically preferred regions. 
Why are masses around 400 GeV required above m^o for a general mixed A^ 12 ? The single most 
important influence causing this characteristic is the aforementioned lack of suppression in the 
thermally averaged cross section, due to N li2 being scalars. 

First, below mw, there are no efficient annihilation channels to deplete N\ except for resonant 
annihilation through the lightest Higgs boson. As this resonance is approached from either above or 
below, the annihilation cross section rises rapidly until the dark matter constraint is satisfied for a 
brief range inm^. This situation is different from other paradigms such as supersymmetry because 
Ni has no direct coupling to the neutral gauge bosons, so there can be no s-channel annihilation 
through Z or 7. Second, once > mw, efficient annihilation occurs to W + W~ . Such significant 
annihilation depresses the relic density below preferred values except for the high-cos 2 region. 
It is this high- cos 2 region where the N\ is primarily an SU(2) singlet, so it cannot annihilate 
effectively to W + W~ . Above mw, the preferred region gradually starts to prefer lower values of 
cos 2, # r/ as m Nl increases. This is because Qdm^ 2 ~ m/(av), so in order to keep fl DM h 2 constant 
for a higher mass, the thermally averaged annihilation cross section must also increase. This trend 
continues until = m h o. 

As was mentioned above, a region below 200 GeV exists with the right relic abundance. This 
region is below the tt and h°h° production thresholds, so annihilation is not extremely efficient 
barring resonances. We display this low-m^ region in greater detail in the right plot of Figure 121 
Below mw , there are no efficient annihilation channels to deplete Ni except for resonant annihilation 
through the lightest Higgs boson. As this resonance is approached from either above or below, the 
annihilation cross section rises rapidly until the dark matter constraint is satisfied for a brief range 
in mjVi- This situation is different from other paradigms such as supersymmetry because N± has 
no direct coupling to the neutral gauge bosons, so there can be no s-channel annihilation through 
Z or 7. This can be seen from the explicit expression for the covariant derivative, Eq. |^ Second, 
once mjvj > mw, efficient annihilation occurs to W + W~ . Such significant annihilation depresses 
the relic density below preferred values except for the high-cos 2 "d^ region. It is this high-cos 2 "d^ 
region where the Ni is primarily an 577(2) singlet, so it cannot annihilate effectively to W + W~. 
Above mw, the preferred region gradually starts to move towards lower values of cos 2 fi^ as m^ 
increases. This is because Q DM h 2 ~ m/(av), so in order to keep VL DM h 2 constant for a higher 
mass, the thermally averaged annihilation cross section must also increase. The trend continues 
until m^x — m^o. Finally, at m^o, the h°h° production threshold is crossed and the relic density 
is uniformly reduced below a relevant level. If one can imagine this parameter space absent of 
the thresholds and resonances, a general trend towards heavier masses at lower cos is apparent 
if m Nl > raw- This continues for the entire parameter space and is due to the strong costi^ 
dependence in the couplings of iVi )2 to the Standard Model W ± . 

Between 200 and 400 GeV in Figure ^ a region is visible without a significant relic density. This 
is due primarily to efficient annihilation into tt and h°h° final states. Preferred regions arise again 
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above 400 GeV and move towards the weak end of the mixed region (low cos 2 ^^) by 800 GeV. 
This region between 400 and 800 GeV is dominated by annihilations to W + W~, so the cosi?^ 
dependence comes directly from mixing angle effects on the gauge couplings of A^ 12 . The main 
additional structure in this region is the slight kink at 500 GeV which is due to the production 
threshold for heavy Higgs final states. The region above 1000 GeV is excluded experimentally 
because annihilation is not efficient enough to reduce the relic density appropriately. It is possible 
that UV completions to this model may provide additional annihilation channels. However, the mass 
of the dark matter particle eventually becomes limited by partial-wave unitarity considerations and 
one has a difficult time avoiding the upper dark matter bound for a weakly coupled thermal relic 
of mass above 2 TeV jH3 EH] • 

How sensitive are these regions of parameter space to the assumptions we have made? The most 
robust regions are the super- weak area between — 80 GeV and — 140 GeV and the sliver 
of parameter space below m h o/2, shown in detail in the right plot in Figure 121 If one assumes that 
the b quark couples proportional to cot (3 instead of tan f3, then the entire region between and 
raw becomes completely depleted of relevant dark matter, but the region between mjVi — 80 GeV 
and m^i — 140 GeV remains unscathed and the cos 2 ^^-independent sliver below remains, 
but drops down to around 50 GeV. Additionally, if one takes a very small value for tan /3, say 
around 0.1, then the preferred band in the region between 400 GeV and 1000 GeV in Figure Q gets 
pushed up to values between 1700 GeV and 2000 GeV, though the existence of the band is robust. 
So, regardless of the assumptions that are made, regions with an acceptable relic density exist below 
m/jo/2 and between 80 and 140 GeV for large values of cos 2 ^^ and also above 900 GeV for any 
value of cos 2 1?^. 

Dark matter direct detection limits depend on the mass of the relic particle. Direct detection 
experiments are designed to best detect particles roughly in the m Nl = 80 — 140 GeV range 4 . In 
this model, a mass in this range requires a high value of cos 2 t?^. To achieve this, it requires a large 
gauge contribution to the (f) mass which translates into heavier super-weak vector bosons. This is 
the limit that precision electroweak observables prefer [7] . A study of direct and indirect detection 
prospects is beyond the scope of this paper, but is the subject of current ongoing work [S|. 

Since gauge contributions dominate the mass matrix (Eq. EJ) for the neutral (r], <fr) system, Eq. 
implies that large cos 2 fi^ happens easily. To investigate this region in more detail, we have focused 
on large cos 2 in the left plot in Figure |21 Taking this plot to represent the 'natural' portion 
of the full parameter space, we can see that over half of the region (the part with red coloring) 
is ruled out on cosmological grounds. Most of the rest of the parameter space has a relic density 
too low to be cosmo logically relevant. The regions with acceptable relic density can be grouped 
into two categories. The first region occurs between 400 GeV and 900 GeV and is due mainly to 
final states W + W~, ti and h°h°. This region does not require any tight relationship between Aq^ 
and another particle. The second category is the low-m^ region below = 200 GeV. If, for 

4 Direct detection possibilities for the super-weakly coupled dark matter candidate are suppressed compared to 
typical supersymmetric candidates for two important reasons. First, there are no intermediate colored states that 
couple the candidate to nuclear matter. Second, there are no couplings between Ni t 2 and the neutral weak vector 
bosons. Thus, nuclear scattering will occur dominantly through s-channel Higgs bosons, which is suppressed by 
Yukawa couplings. 
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Figure 2: Relic Density Contours for LH Models. This Figure is the same as Figure ^ but 
magnified to show the finer detail in the low mass region. The left plot also focuses in on the region 
with large cos 2 ^^. 



example, we take cos 2 d v< f> = 0.95 and review Figure El we find that the dark matter candidate must 
take very specific mass values: mjvi = "m>h /^ ± 5, mw, n^h , or 400 ± 50 GeV. Again, 5 is a small 
number around 3 GeV. While it is possible for cos 2 $^ to be large without being around 0.95, the 
predictive power of such a scenario in terms of dark matter is intriguing. 

In summary, we find that LH models with a general mixed Ni state do not provide enough relic 
density to explain dark matter unless either additional Q plaquettes are included to lift m^, or 
m Nl is chosen to have very specific values, namely mw or m h o/2 ± 3 GeV. However, we do find 
that the most likely scenario, that of a super-weakly interacting Aq with ~ 100 GeV, is able 
to provide for the correct relic density, although the preferred region is quite narrow. 

3.1 Other Possible Effects 

At this point it is natural to inquire about other contributions to the relic density calculations 
that might change our results. As we will describe, additional effects will likely only further reduce 
the available parameter space. Possible contributions fall into one of three categories: additional 
annihilation channels, coannihilation channels, and entropy generation after Aq freeze-out. We will 
discuss each of these in turn. 

There are almost certainly other ways that two Aq particles can annihilate. These new channels 
may come from other loop-suppressed operators, or possibly from operators that appear in other 
specific models. Any time new annihilation channels are allowed, one expects the thermally averaged 
cross section to be increased. This will further suppress the relic density of Aq. We have already 
seen this in Figures Q and El when the resonant annihilation and particle production thresholds were 
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crossed. 



As mentioned previously, it is also possible for the mass of <pf (the lighter of 0+ and 0+) to 
be degenerate with the mass of Ni to within 5-10%. Such a scenario can occur when the mass 
contributions coming from the f2-plaquettes dominate the mass matrices. While this situation is 
unlikely for light Aqs, it is not unnatural for heavy Aqs, whose masses can be principally determined 
by fi-plaquettes. When this degeneracy happens, the two particles freeze out of equilibrium at 
about the same time, and the interactions of <f)f can significantly affect the resulting relic density 
of Ni HE] • We can make the normal assumption that <j)f decays into Aq on a time scale much 
shorter than the Hubble expansion time scale. Thus, instead of treating the evolution of the number 
density rii of each species separately, we can evolve the total number density n = n Nl + n^+. To do 
this, we need to approximate the thermally averaged cross section as 



where the eq superscript denotes values at equilibrium and i and j denote particle species to be 
summed over (Aq and (pf for this example). From this formula we can see that just adding a 
coannihilator with equal annihilation cross sections to the original particle has little effect on the 
relic density. In fact, this is the case for our analysis already. We have included the fact that 
t] x and 7]y (as well as <p x and <p y ) are very nearly degenerate in mass and can coannihilate. In 
order for coannihilation to have a significant effect on the relic density, it is necessary for some of 
the coannihilation channels to dominate over the normal annihilation channels. This happens in 
mSUGRA when the lightest neutralino can coannihilate with the lightest stau |27| 128] or stop |29j . 
It is also possible to have coannihilation with another neutralino ^T] both in mSUGRA in the 'focus 
point' region (SOj and also in rSUGRA [HI] and string derived models [32] . In LH models, the net 
effect is expected to be qualitatively similar to supersymmetric models. Including coannihilation 
effects, like including additional annihilation channels, will reduce the relic density. 

It is also possible that significant entropy can be generated between the freeze-out temperature 
of Aq and the present temperature. When the relic density is calculated, the quantity that is evolved 
is the number density per comoving volume. This is computed by dividing the number density by 
the entropy density: — n^/s. In making this redefinition, it is assumed that entropy is neither 
added to nor removed from the system under study. However, such an assumption about the entropy 
need not hold. If entropy is added to or taken away from the universe, then the number density per 
comoving volume will not obey the simple standard evolution. In terms of an existing relic density, 
the net effect will be to multiply a particle's relic density by a factor of roughly S u/ S new . S id 
is the entropy per comoving volume before the entropy change occurred and S new is the entropy 
per comoving volume after the change in entropy occurred. The end result is again a reduction 
in the relic density for the case of entropy generation 5 . In some sense this generic reduction in 
relic density is good because it also reduces the amount of parameter space that is excluded on 
cosmological grounds (where fldmh 2 > 0.129). However, it also results in a reduction in the amount 

It is possible, however, to have both entropy production and relic density increase through the decay of a heavy 
modulus |531I53] . 




(29) 
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of cosmologically preferred parameter space (where 0.094 < VL dm h 2 < 0.129). Thus, only the regions 
with Qdmh 2 < 0.094 are enlarged. Such regions do not supply enough relic density from the N\ t 2 
particles, so one might have to appeal to dark matter candidates outside of the LH model to meet 
current dark matter density measurements. 

4 Conclusion 

We have investigated the possibility of having a dark matter candidate arise from a 'little Higgs' 
model. We found regions of parameter space which satisfy the relic density requirement 
0.094 < Qdmh 2 ^ 0.129. These regions require a dark matter candidate that either has a mass 
between roughly m/jo/2 and m^o or has a mass above 400 GeV. The low-mass parameter space 
additionally requires the candidate to have either 'super-weak' coupling or very specific masses of 
myfijl ± 3 GeV or m\y. This is somewhat different from the standard dark matter scenario in 
theories involving broken supersymmetry and stems from the LH dark matter candidate being a 
scalar and also not having any direct couplings to the neutral vector bosons. We discussed additional 
possible contributions to the relic density calculation of the dark matter candidate. We found that 
the additional contributions tended to result in an even lower relic abundance. 

The standard estimate for dark matter gives that weak couplings and masses O(10 GeV) — 
(9(1 TeV) result in the right order of magnitude for the dark matter abundance. This numerology 
indicates that dark matter might be related to electroweak symmetry breaking or the physics 
stabilizing the weak scale. In this note we have shown that it seems possible to accommodate 
sufficient dark matter for theory space models with exact discrete symmetries. If this had not 
been the case either by over-depletion or by a breaking of the discrete symmetry, then there would 
still be the usual dark matter problem. Through standard calculations using partial-wave unitarity 
it is possible to put an upper limit of 550 TeV on the mass of the dark matter particle ! 35 . By 
using a more refined calculation assuming weak coupling, the limit can be brought down to about 
1.8 TeV (US]. If there is insufficient dark matter from such theory space models and if the dark 
matter is related to the weak scale, then it must come from physics in the ultraviolet and the 
particles must have masses roughly O(f) or less if they are weakly coupled. However, it should also 
be noted that it is possible that dark matter can originate from axions or also a sector unconnected 
to the Standard Model except through gravity. 

In conclusion, we have found viable dark matter in a prototypical LH model. Requiring the 
correct relic density forces the LH model to one of three regions. First, iV lj2 can have masses of 
m/jo/2 ± 3 GeV or m\y, resulting in resonant and threshold annihilation, respectively. Second, it is 
possible to have a 'super-weak' AT 1)2 with cos 2 1?^ ~ 1 (where the gauge contributions dominate in 
the rj-(p mass matrix). Lastly, the LH model can have large Q plaquettes to raise mN l2 - However, 
both the extreme 'super-weak' regime and very large contributions from Q plaquettes appear to 
over-close the universe. 
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